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Abstract 

A coupled channel analysis of the centrally produced K^K~ and tt^tt^ final states has been 
performed in pp collisions at an incident beam momentum of 450 GeV/c. The pole positions 



^ 

Q I and branching ratios to tttt and KK of the /o(980), /o(1370), /o(1500) and /o(1710) have been 

P_i' determined. A systematic study of the production properties of all the resonances observed in 

D ' the TT+TT" and K^K~ channels has been performed. 
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Recently the WA102 collaboration has published the results of partial wave analyses of 
the centrally produced K^K~ , K^K^ vr+vr" and vr^vr'^ channels. A striking 
feature of these analyses was the result that the /j(1710) has J = (we shall refer to it as the 
/o(1710) hereafter). In these papers the S-wave from each channel was fitted independently 
using interfering Breit-Wigners and a background. In this present paper we will first show 
how the resulting parameters change if a different method of fitting is used, namely, a T- 
Matrix analysis and a K-Matrix analysis using the methods described in ref. Next we will 
perform a coupled channel fit to the K^K~ and tt^tt' final states in order to determine the 
pole positions and branching ratios of the observed mesons. Finally we will present information 
on the production kinematics of these resonances. 

In our previous publication a fit has been performed to the 7r"'"7r~ S-wave using a coherent 
sum of relativistic Breit-Wigner functions and a background of the form: 

n=l 

where the background has been parameterised as 



7V7V 



where a„ and 9n are the amplitude and the phase of the n-th resonance respectively, a, /3, 
7 and 6 are real parameters, BW{Mt^t^) is the relativistic Breit-Wigner function for a spin 
zero resonance. In order to describe the centrally produced ti^tt~ mass spectrum the function 
|A(M7r7r)P has been multiplied by the kinematical factor {Mt,.„ — Am^Y^"^ / M^^ 0. The resulting 
function is then convoluted with a Gaussian to account for the experimental mass resolution. 

In this present paper we use the Flatte formula to describe the /o(980), this is referred 
to as Method I. For the tt+tt" channel the Breit-Wigner has the form: 



mo — — ■jmo(r^ - 
and in the K^K~ channel the Breit-Wigner has the form: 
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BW{Mkk) 



where Fj is absorbed into the intensity of the resonance. F^^ and F^ describe the partial widths 
of the resonance to decay to vrvr and KK and are given by 



K 



where g-^ and gx are the squares of the coupling constants of the resonance to the vrvr and KK 
systems. The resulting fit is shown in fig. ^) for the entire mass spectrum and in fig. ^d) for 
masses above 1 GeV. The sheet II pole positions [0 for the resonances are 

/o(980) M = ( 983 ± 8) -z ( 58 ± 11) MeV 

/o(1370) M= (1306 ±18) -i(lll±23) MeV 

/o(1500) M = (1502 ± 12) -i ( 65 ± 12) MeV 

/o(1710) M= (1748 ±22) ( 73 ± 22) MeV 
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These parameters are consistent with the PDG values for these resonances. 



To test the sensitivity of these resuhs on the fitting method used we have also performed a 
fit to the TT+TT" mass spectrum using the T-Matrix parameterisation of Zou and Bugg |@. The 
invariant amplitude for tt^tt^ central production can be expressed as 

A = ai(s)Tii + a2{s)T2i (1) 

where Tn and T21 are the invariant amplitudes for elastic nn nir and KK — > nn scattering 
and are parameterised by 

2ipi Ml- s- i{pigi + p2^2) 
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(3) 



21 — 

Mr- s- i{pigi + p2g2) 

where pi = (1 — 'im^./sY^'^ and p2 = (1 — 4m|^/s)^/^ are phase space factors and s is the 
invariant mass squared of the tt+tt^ channel. The background term is presumed to be coupled 
only to the vrvr channel and has the form 



e 



2*0 _ ]^ I 



2ipi A{s) — ipi 

which satisfies the unitarity condition. A{s) is an arbitrary real function which has been taken 
to be of the form 

... 1 + ais + a2S^ 

6i(s-m2/2)+62s2 

The real functions ai{s) in equation (|l|) describe the coupling of the initial state to the 
channel i. These functions are approximated by the power expression 0: 

„=0 

where the factor Amj^ is introduced as a convenient scaling. It has been found that n = 3 is 
sufficient to describe the S-wave distribution. 

In order to describe the centrally produced tt+tt^ mass spectrum, the function |A(M7r7r)P has 
been multiplied by the kinematical factor (M^^r — 4m^)^/^/M^^ and the resulting function is 
then convoluted with a Gaussian to account for the experimental mass resolution. The resulting 
fit is shown in fig. |T]c) for the entire mass spectrum and in fig. |I|d) for masses above 1 GeV. 
As can be seen the fit does not describe well the region above 1.0 GeV. The sheet II pole 
corresponding to the /o(980) is 

M = ( 993 ± 8) -i ( 38 ± 9) MeV 

There are two poles from the background term with 

Ml = ( 388 ± 55) -i (223 ± 28) MeV 
M2 = (1541 ± 32) -i (143 ± 21) MeV 
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The first pole may be associated with the low mass S-wave. The NA12/2 collaboration |^ 
have previously shown that this region may be interpreted as being due to the a particle. The 
second pole would appear to be in the region of the /o(1500) but the width is very broad. This 
is due to the fact that the fit is not able to properly describe the region around 1.3 GeV. 

Adding one more term to equations (H) and (§) to describe the 1300 MeV mass region 
improves the fit considerably but still does not describe the region around 1700 MeV. In order 
to produce a satisfactory fit two terms have to be added to equations (H) and @ to account 
for the /o(1370) and /o(1710) which results in the fit shown in fig. |I|e) for masses above 1 GeV. 
The new sheet II pole positions are 



/o(980) 


M 


= ( 992 


± 


6) 


-i( 52 ± 9) 


MeV 


/o(1370) 


M 


= (1310 


± 


30) 


-i(134 ± 23) 


MeV 


/o(1500) 


M 


= (1497 


± 


17) 


-il 82 ± 21) 


MeV 


/o(1710) 


M 


= (1752 


± 


15) 


-il 53 ± 12) 


MeV 



These parameters are consistent with the values from the fit using interfering Breit-Wigners 
and with the PDG |^ values for these resonances. 

An alternative parameterisation is to use the K-Matrix formalism. In this case the Lorentz 
invariant T-Matrix is expressed as 



f = k{i-ipk) 
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where for the case of vrvr and KK final states p is a 2 dimensional diagonal matrix and 7^ is a 
real symmetric 2 dimensional matrix of the form 

- s Sb- s 

In the K-Matrix formalism the background is assumed to couple to both the vrvr and KK chan- 
nels. In order to describe the centrally produced tt+tt" mass spectrum, the function |y4(Mjr^)p 
has been multiplied by the kinematical factor (Mtttt— 4m^)^/^/M^^ |^ and the resulting function 
is then convoluted with a Gaussian to account for the experimental mass resolution. 

Two coupled channel resonances are found from the fit shown in fig. |l]f ) for the entire mass 
spectrum and in fig. |l^) for masses above 1 GeV with their sheet II T-Matrix poles at 

Ml = ( 988 ± 18) -i ( 39 ± 12) MeV 
Ms = (1526 ± 22) -i(191 ± 53) MeV 



As in the case of the T-Matrix analysis the fit fails in the 1.3 GeV region. Adding one additional 
pole improves the fit. However, in order to to obtain a satisfactory fit it is found necessary 
to include two extra poles. The fit shown is in fig. |I]h) for masses above 1 GeV and results in 
sheet II pole positions of 

/o(980) M = ( 982 ± 9) -i ( 38 ± 16) MeV 

/o(1370) M = (1290 ± 30) -i(104 ± 25) MeV 

/o(1500) M = (1510 ± 10) -i l 56 ± 15) MeV 

/o(1710) M = (1709 ± 15) -i l 75 ± 18) MeV 
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These parameters are consistent with the values from the two previous fits. 



Finally, in order to perform a coupled channel fit to the vr+vr" and K^K'^ final states a 
correct normalisation of the two data sets has to be performed. The fit has been modified 
to take into account the relative differences in geometrical acceptance, event reconstruction 
and event selection. The fit also includes corrections for the unseen decay modes so that 
the branching ratio vrvr to KK can be calculated. In addition to the resonances discussed 
above, the ao(980) can also contribute to the KK S-wave mass spectrum. The contribution 
from ao(980) — > K^K~ has been calculated from the observed decay ao(980) rjn and the 
measured branching ratio of the ao(980) [T^. The calculated contribution (500 ± 120 events) 
has been included in the fit as a histogram. There is no evidence for any ao(1450) contribution 
in the rjn channel and hence it has not been included in the fit to the KK S-wave. 

The coupled channel fit has been performed using the three methods described above. The 
following pole positions and branching ratios quoted for the resonances are a mean from the 
three methods. The statistical error is the largest error from the three fits and the systematic 
error quoted represents the spread of the values from the three methods. The results of the 
combined fit are shown in fig. |^. The sheet II pole positions are 

/o(980) M = ( 987 ± 6± 6) -i ( 48 ± 12 ± 8) MeV 

/o(1370) M = (1312 ± 25 ± 10) -i(109 ± 22 ± 15) MeV 

/o(1500) M = (1502 ± 12 ± 10) -i( 49 ± 9 ± 8) MeV 

/o(1710) M = (1727 ± 12 ± 11) -i( 63 ± 8 ± 9) MeV 

These parameters are consistent with the PDG @] values for these resonances. For the /o(980) 
the couplings were determined to be = 0.19 ± 0.03 ± 0.04 and gx =0.40 ± 0.04 ± 0.04. 

The branching ratios for the /o(1370), /o(1500) and /o(1710) have been calculated to be: 

/o(1370) ^ KK 



/o(1370) ^ vrvr 

/o(1500) -^KK _ 
/o(1500) ^ vrvr 

/o(1710) ^ KK 



/o(1710) 



0.46 ±0.15 ±0.11 

0.33 ±0.03 ±0.07 
= 5.0 ±0.6 ±0.9 



TTTT 



These values are to be compared with the PDG values of 1.35 ± 0.68 for the /o(1370) 
and 0.19 ± 0.07 for the /o(1500), which comes from the Crystal Barrel experiment [|l^. The 
value for the /o(1710) is consistent with the value of 2.56 ± 0.9 which comes from the WA76 
experiment which assumed J = 2 for the /j(1710) ||12|| . 

In our previous publications a study has been performed of resonance production rate as 
a function of the difference in the transverse momentum vectors {(IPt) between the particles 
exchanged from the fast and slow vertices |l^, [l^. It has been observed ||15|| that all the 



undisputed qq states (i.e. rj, t]', /i(1285) etc.) are suppressed as (IPt goes to zero, whereas 
the glueball candidates /o(1500) and /2(1950) survive. In order to calculate the contribution of 
each resonance as a function of dPT, the partial waves have been fitted in three (IPt intervals 
with the parameters of the resonances fixed to those obtained from the fits to the total data 
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using method I. As an example of how the mass spectra change as a function of (IPt, fig- ^ 
shows the tt+tt^ S-wave and D-wave in three (IPt intervals. 

The S-wave clearly shows that the region around 1.3 GeV is more enhanced at large (IPt 
in contrast to the low mass part of the spectrum. This is not due to feed through from the 
D-wave which has been estimated to be less than 3 % in this mass region. In the D-wave the 
/2(1270) is suppressed at small (IPt in contrast to the behaviour of the low mass part of the 
D-wave. 

Table |l] gives the percentage of each resonance in three (IPt intervals together with the ratio 
of the number of events for (IPt < 0.2 GeV to the number of events for (IPt > 0.5 GeV for 
each resonance considered. As can be seen from table |l|, the p°(770), /2(1270) and /2(1525) are 
suppressed at small dPT in contrast to the /o(980), /o(1500) and /o(1710). 

The azimuthal angle (0) is defined as the angle between the pt vectors of the two protons. 
In order to determine the dependence for the resonances observed, the partial waves have 
been fitted in 30 degree bins of with the parameters of the resonances fixed to those obtained 
from the fits to the total data. The fraction of each resonance as a function of is plotted in 
fig. ^ The dependences are clearly not flat and considerable variation is observed between 
the different resonances. 

In order to determine the four momentum transfer dependence (t) of the resonances observed 
in the 7i~^7c~ and K^K~ channels the partial waves have been fitted in 0.1 GeV^ bins of t with 
the parameters of the resonances fixed to those obtained from the fits to the total data. Fig. ^ 
shows the four momentum transfer from one of the proton vertices for these resonances. The 
distributions for the /o(980), /o(1370), /o(1500), /o(1710) and p(770) have been fitted with a 
single exponential of the form exp{—h\t\) and the values of h found are given in table ^ The 
distributions for the /2(1270) and /2(1525) cannot be fitted with a single exponential. Instead 
they have been fitted to the form 

dt 

The parameters resulting from the fit are given in table ^. 



In a previous publication by the WA76 collaboration |T6| the ratios of the production cross 
sections for the p(770), /o(980) and /2(1270) were calculated at = 12.7 and 23.8 GeV. 
However, the experiment at 300 GeV (^i = 23.8 GeV) was only sensitive to angles less 
than 90 degrees and the acceptance program that had been used assumed a flat distribution. 
Hence the cross sections at 300 GeV were underestimated for the p(770) and /2(1270) and 
overestimated for the /o(980). 

After correcting for geometrical acceptances, detector efficiencies, losses due to cuts, and 
unseen decay modes, the ratios of the cross-sections for the p(770), /o(980), /2(1270) and in 
addition, the /o(1500) at ^/s = 29.1 and 12.7 GeV are given in table The cross sections for 
the /o(980), /2(1270) and /o(1500) are compatible with being independent of energy which is 
consistent with them being produced via double Pomeron exchange Jl^ . 



In summary, a coupled channel fit has been performed to the centrally produced tt+tt and 
K^K~ mass spectra. The pole positions and branching ratios of the /o(980), /o(1370), /o(1500) 
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and /o(1710) have been determined using three different methods which give consistent results. 
An analysis of the dP^ dependence of the resonances observed shows that the undisputed qq 
mesons are suppressed at small dP^ in contrast to the enigmatic /o(980), /o(1500) and /o(1710). 
Considerable variation is observed in the 4> distributions of the produced mesons. 
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Table 1: Production of the resonances as a function of (IPt expressed as a percentage of tlieir 
total contribution and the ratio (R) of events produced at dPx < 0.2 GeV to the events produced 
at dPx > 0.5 GeV. 





dPrr<Q 2 GeV 


2<dPT<0 5 GeV 


dP^>Q 5 GeV 


p dPT<0.2GeV 
~ dPT>0-5GeV 




oo 1 _|_ 9 9 


ccfl « _|_ 9 7 


9« 9 _|_ 9 7 


n CO _|_ f) 1 9 




1 Q 1 -1-4 


090190 


4Q Q _|_ 9 Q 


U.OU ZIZ u.uo 






C4 1 -1- c: n 


99 _|_ f) 

ZiZi.O ZIZ o.u 


1 nc _|_ fi 10 
x.uu ZIZ yj . ±0 






4 5 Q _)_ 1 n 

THj . _l_ -L . W 


27 s + 1 1 

^1.0 1 J- . J- 


n 05 + n nfi 

W.tytJ _l_ \J.\J\J 


p(770) 


6.4 ± 2.5 


41.1 ± 4.0 


52.5 ± 3.0 


0.12 ± 0.05 


/2(1270) 


7.6 ± 1.2 


29.5 ± 0.6 


62.9 ± 1.7 


0.12 ± 0.02 


/^(1525) 


4.3 ± 3.0 


35.7 ± 3.0 


60.0 ± 4.0 


0.07 ± 0.04 


/2(2150) 


2.6 ± 3.0 


53.6 ± 3.6 


43.6 ± 3.2 


0.06 ± 0.08 
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Table 2: The slope parameter b from a single exponential fit to the \t\ distributions. 



/o(980) 


/o(1370) 


/o(1500) 


/o(1710) 


P(770) 


b/GeV-2 5.3 ± 0.2 


6.7 ±0.5 


5.2 ±0.5 


6.5 ± 0.8 


5.8 ±0.1 



Table 3: The slope parameters from a fit to the \t\ distributions of the form ^ — ae °'^^+(3te 





a 


bi 




&2 






GeV-2 




GeV-2 


/2(1270) 


0.25 ± 0.04 


8.3 ±3.0 


5.3 ±1.0 


8.7 ±0.4 


/^(1525) 


0.4 ±0.3 


8.3 ±5.0 


5.3 ±4.0 


7.7 ±3.3 



Table 4: The ratio of the cross sections at ^/s — 29.1 and 12.7 GeV. 



p(770) 


/o(980) 


/2(1270) 


/o(1500) 


fJ^Sfi 0.36 ±0.05 


1.28 ±0.21 


0.98 ±0.13 


1.07 ±0.14 
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Figures 

Figure 1: The tt'^tt' mass spectrum with a) and b) fit using interfering Breit-Wigners. Fit 

using the T-Matrix parameterisation c) and d) the original fit and e) including two extra poles. 
Fit using the K-Matrix parameterisation f ) and g) the original fit and h) including another two 
poles. 



Figure 2: A coupled channel fit to the n'^n and K'^K S-wave distributions. 



Figure 3: The tt+tt S-wave and D-wave in three dPr intervals. 



Figure 4: The distributions of the resonances observed in the tt+tt and K^K channels. 



Figure 5: The four momentum transfer squared (|t|) from one of the proton vertices for the 
resonances observed in the 7r"'"7r~ and K^K~ channels. 
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